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1 Introduction 



The PIBETA project at the Paul Scherrer Institute (PSI) is a program of mea- 
surements which aims to make a precise determination of the tt + — > 7r°e + z/ e 
decay rate, initially with an uncertainty of ~ 0.5% [1,2]. The small branching 
ratio (~ 1 x 10~ 8 ) for the pion beta (ir/3) decay and the required high measure- 
ment precision impose stringent requirements on the experimental apparatus. 

The experimental signature of a 7r + — > 7r°e + z/ e event is determined by the 
prompt decay 7r° — > 77. The detector must be able to handle high event rates 
and cover a large solid angle with high efficiency for n° detection. Efficient 
hardware suppression of background events requires good energy and timing 
resolution. At the same time the system must operate with low systematic 
uncertainties and be subject to accurate calibration. 

We have chosen to detect nf3 decays at rest, and to use n + — > e + u e (7r e2 ) 
decays for normalization of the absolute branching ratio. Consequently, our 
apparatus has the following main components [3]: 

(1) a forward beam counter (BC), an active degrader (AD) and a segmented 
active target (AT) to stop the beam pions; 

(2) two concentric cylindrical multi-wire proportional chambers (MWPCi and 
MWPC2) surrounding the active target for charged particle tracking; 

(3) a segmented fast veto scintillation counter (PV) surrounding the MWPC's, 
for the particle identification; 

(4) a high resolution, 240-module, pure Csl segmented fast shower calorimeter 
surrounding the active target and tracking detectors in a near-spherical 
geometry; 

(5) cosmic ray veto counters (CV) around the entire apparatus. 

The 7i/3 branching ratio B^p is evaluated from the following expression: 

D f^ £ 2 ' Bw e 2 A ne2 e 7re2 . . 

B^p - —b j~r a { > 



where B w p, B We2 and B 1 ^^ ni are the branching ratios for the 7r(3, 7r + — * e + v e 
and 7T° — > 77 processes, respectively, N^p and N ne2 represent the detected 
numbers of good nf3 and 7r + — > e + z/ e events, A n p and A ne2 are the detector 
acceptances, e^p and e ne2 are the detector efficiencies, and f ne2 is the prescaling 
factor for the 7r + — > e + u e trigger. 

Due to the similarities between the two classes of events, the acceptances A n p 
and A We2 are nearly equal. At the same time, both are functions of 
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(1) the geometrical solid angle covered by the calorimeter; 

(2) geometry of a beam stopping distribution; 

(3) energy line-shapes of photon- and positron-induced showers in the calorim- 
eter; 

(4) software cuts imposed on the measured event variables. 

The main source of uncertainty of geometrical nature in the absolute accep- 
tance ratio A^ e2 /A^p has to do with the position and spatial spreading of the 
pion beam stopping distribution. Due to the cylindrical symmetry of our detec- 
tor system the uncertainties can be broken down into lateral (x-y coordinates 
perpendicular to the beam direction) and axial uncertainties (z coordinate in 
the direction of the beam axis). Of the two, the lateral uncertainties produce 
a larger effect on the acceptances than the axial uncertainties. Also, since the 
n P events require symmetrical coverage for the two photons emitted in the n° 
decay (nearly at rest), A n p is affected more than A ne2 . 

In this report we analyze the PIBETA data we have collected in order to mon- 
itor the centroid and shape of the pion beam stopping distribution. Our goal 
is to examine how the uncertainty in the deduced beam geometry affects the 
A We2 /A 7rj 3 ratio. In Section 2 we specify the design and layout of the experi- 
mental beam line. Sections 3 and 4 address the related issue of geometrical 
uncertainties associated with the assembly of the modular shower calorimeter. 
In Section 5 we describe the PIBETA detector's central region, most impor- 
tantly the segmented active target. Section 6 lays out the method and specifies 
the accuracy with which the lateral distribution of stopping pions (and con- 
taminating positrons) can be deduced from their respective calibrated pulse 
height spectra in the active target segments. Section 7 introduces an alter- 
native way of determining the beam position and shape that relies on the 
back-tracking of charged particles using MWPC hits. In Section 8 the beam 
stopping distribution is deduced from the reconstructed vertices of the Dalitz 
decay events tt° — > 7e + e~ which originate from the pion beta decays at rest. 
Influence of the different beam distributions on (i) the absolute acceptance of 
the detector, and (ii) measured n° — > 77 opening angle distributions, is calcu- 
lated in a Monte Carlo (MC) detector simulation and presented in Section 9. 
Finally, the mutual consistency of the beam profile reconstructions by these 
four independent methods is summarized in the Conclusion. 



2 PIBETA detector beam line 



The measurements are performed in the nEl channel at PSI [4]. The beam 
line is operated in the high intensity, low momentum resolution mode. A 114 
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MeV/c 7r + beam tune is developed with momentum spread Ap/p < 1.2 % and 
maximum nominal 7r + beam intensity I n of ~ 1.2 • 10 6 7r + 's per second. 

The spatial spread of the n + beam is restricted by a lead collimator with a 
7mm pin-hole located 3.985 m upstream of the detector center. The beam 
particles are first registered in a 3 mm thick plastic scintillator placed directly 
behind the collimator. The pions are subsequently slowed in a 40 mm long ac- 
tive plastic degrader and stopped in an active plastic target that is positioned 
in the center of the PIBETA detector. The beam line properties have been 
computed with TRANSPORT [5] and TURTLE [6] codes. The TURTLE mo- 
mentum spectrum of 7r + 's incident on the front face of the degrader counter 
is shown in Fig. 1. 

This specific choice of the incident beam momentum allows for the optimum 
pion/muon/positron separation by the time-of-flight (TOF) method. The ab- 
solute TOF's measured between the beam counter and the degrader /target 
detectors for the beam positrons, muons, and pions are 12.7 ns, 17.4 ns, and 
20.2 ns, respectively. 

The beam particle triggers are defined by a four-fold coincidence between 
(1) forward beam counter BC, (2) active degrader AD, (3) active target AT, 
and (4) the accelerator rf signal. The energy depositions of minimum-ionizing 
positrons in BC, AD, and AT counters are 0.6 MeV, 7.2 MeV, and 9.0 MeV, 
respectively. The corresponding energy deposition values for 114MeV/c pions 
(with kinetic energy 40 MeV) are 0.7 MeV, 13.0 MeV, and 28.0 MeV, respec- 
tively. By appropriately adjusting the discriminator thresholds and time delays 
of the logic signals that make up the quad coincidence, the "7r-in-beam" and 
"e-in-beam" triggers, which preferentially select beam pions and positrons, re- 
spectively, are set up. A digital oscilloscope snap-shot of four trigger-defining 
signals at the input of the "7r-in-beam" coincidence unit is shown in Fig. 2. 

The e + and fj, + beam contaminations in the "7r-in-beam" trigger measured in 
the TOF spectra are small, ~0.4% and ~0.2%, respectively. 



3 Design and Assembly of the Csl calorimeter sphere 



The design and assembly of the Csl calorimeter directly influence the precision 
with which the geometry of the beam and the decay particle tracks can be 
determined. 

The spherically shaped pure Csl calorimeter is mounted inside a forged spheri- 
cal steel shell which itself is supported by a steel frame. The shell has two large 
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axial openings of 550 mm diameter for beam entrance and exit and 220 smaller 
holes distributed over its surface and aligned with the axes of the individual 
Csl crystals. The smaller holes were made for access to the Csl photomulti- 
plier tubes (PMT's). The beam openings were used to hold cast iron plugs 
with truncated pyramid shapes inside the shell to provide contact surfaces for 
20 outlying shower veto crystals, respectively. These two plug plates press the 
stacked crystals together and hold them in place. The plugs have a cylindri- 
cal bore-hole (diameter 270 mm) for the beam which also allows access to the 
interior of the Csl sphere, e.g. the wire chambers, the target counters, etc. 

Assembly of the calorimeter, weighing approximately 1,500 kg, proceeded as 
follows: 

The empty spherical shell with its axis turned vertical, supported by its steel 
frame, was equipped with a single plug insert from below. The beam hole of 
this plug was used to hold a micrometer gauge with a touch-ball in place. 
This device allowed the measurement of a crystals' front surface distance to 
the center of the sphere with 100 /im accuracy. Afterwards, the whole frame 
was turned with the plug at the bottom. Through the other beam hole, on 
top now and not yet closed with its plug, all crystals already equipped with 
their PMT's were moved into the interior of the spherical shell. First came 
10 pieces of trapezohedron-shaped veto crystals which were positioned on the 
flat surfaces of the upstanding pentagonal pyramid of the bottom plug. The 
correct distance to the center of the sphere was adjusted using threaded bolts 
mounted on the plug and supporting the crystals' rear end. 

After all 10 bottom veto crystals were in place, the next crystal layer was 
stacked on top of the veto layer. Threaded tubes with bore-holes large enough 
to house a PMT and its mu-metal shielding were installed through the smaller 
holes in the shell, which are aligned with the crystals' axes. The front end 
of those tubes pressed against the rear surface of the crystals and thereby 
allowed for a precise adjustment of their distances to the detector center point, 
as measured with the touch-ball measurement device. Then the high voltage 
divider base was put through the tube and connected with the PMT. Bases for 
the veto crystals had to be mounted somewhat differently as there are no holes 
available for them on the shell due to the close proximity of the axial plugs. 
The closest hole, later employed for a crystal on the next layer, was used to 
move the veto crystals base into the sphere, then sidewards until connection 
the the veto's PMT was possible. Small supporting structures were mounted 
to press the base gently against the PMT to prevent loosening due to gravity's 
pull for vertically mounted crystals. 

When the sphere was half filled with crystals the procedure had to be changed 
somewhat as the upper surface of completed layers became steeper and steeper. 
In order to prevent the crystals from sliding down towards the sphere's center, 
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the measurement device at the center was replaced by a pin-holder with a pin 
for each crystal still to be inserted and already adjusted in length to give the 
correct distance of 260 mm to the center. 

After the last layer of 10 veto crystals was filled in, the second pentagonally 
shaped plug was mounted on top and pressed against the crystal ball under- 
neath. Calculations made prior to the assembly told us that pressing the plug 
by 10 mm against the crystal ball would be needed in order to pack all crystals 
tightly. This prediction turned out to be exactly fulfilled, giving us confidence 
in the correct positioning of all modules. The completed setup was turned 
around several times and watched for stability before the pin-holder at the 
center was disassembled through the beam entry/exit holes. 

The detector geometry and location has been measured by the PSI surveying 
group using a theodolite. The upstream and downstream detector side are lev- 
elled, with the center 1506.0 ±0.3 mm above the floor, the value corresponding 
to the nominal beam height in the 7rEl channel. The measurement of the de- 
tector center was reproducible and the alignment of the target axis along the 
detector longitudinal axis is confirmed with 0.3 mm uncertainty. 



4 Alignment of the MWPC's and the Csl calorimeter 



Two out of the four reconstruction methods described in this paper depend 
on the information provided by two independent detectors, namely the read- 
ings of the MWPC's and the signals from the Csl calorimeter. Great care was 
taken in the proper assembly and positioning of these two detectors as already 
described in the previous section. During data acquisition, however, wire cham- 
bers required several maintenance sessions and they were once removed from 
the entire assembly and then reinserted. This could have introduced different 
displacements between the wire chamber centers and the center of the Csl 
sphere. Such a displacement could significantly affect both the reconstruction 
of the beam's stopping distribution as well as the measured n° — > 77 open- 
ing angle distribution. The following section describes how we ascertained the 
relative offsets between the centers of the calorimeter and the wire chamber 
coordinate systems. 

On a single event basis it is not possible to find the exact initial point of inter- 
action between a positron/photon and the face of the calorimeter relying only 
on the calorimeter signals. We can, however, determine whether a particular 
Csl detector was close to a given track. Assuming that in such most of 

the energy was deposited in the intersecting crystal, we can adopt the follow- 
ing procedure. We accept events with a single charged particle (e.g. \i — > ev^U^ 
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positron) and reconstruct the MWPC track and its intercept with the calorim- 
eter. We then scan all the energies deposited in the calorimeter modules and 
find the maximum value. Calculated values of the intercept coordinates cor- 
responding to the crystal with most energy deposited can be histogrammed. 
The outcome of this procedure is a set of individual crystal maps, e.g., two- 
dimensional plots (polar angle 9 versus azimuthal angle 0) of a given crystal 
"illumination" as seen by the wire chambers. Knowing the geometry of the 
crystal we can predict how polar angle 6 and azimuthal angle projected 
histograms should look and then compare them to the experimental data. 

Since the inner radius of the calorimeter sphere is fixed at 260 mm we have 
effectively reduced our problem to a two-dimensional one. Namely, everything 
depends only on the spherical coordinates 9 and <f> of the crystals. We know 
that Michel events are uniformly distributed in solid angle. Therefore, if in 9- 
(f> space our crystals were rectangular, then the probability density of hitting 
one particular crystal (for, say, the 9 projection) would be dA/d9 (where dA 
is infinitesimal area of the 0-(f> projection) and therefore would be a constant. 
In 9-(j) space our crystals do not have rectangular shapes and thus probability 
density is not a constant. 

Knowledge of the exact geometry of the crystals allows us to calculate the the- 
oretical dA/d9 analytically. The complete derivation is presented in Ref. [7]. 
Our 9-(j) plots are additionally smeared by the chamber resolution. Mathe- 
matically, it means that the dA/d9 function is convoluted by some resolution 
function. In Sec. 7.2 we show that appropriate MWPC response is a Gaussian 
convolution and therefore the trial function is 



[9 — my 



f( \ ( dA ^ r 



a 2 



d9, (2) 



where T, m and a are the parameters of the Gaussian distribution, determined 
in a fit. 

The exact position of the crystal enters the fitting function through the ana- 
lytical expression for dA/d9 as well as via the limits of integration. 

Fitting this function to the measured 9 projection of the crystal map gives the 
position of the crystal vertices as seen by the chambers. A simple geometric 
transformation allows us to calculate the longitudinal (z) displacement of the 
chambers with respect to the Csl sphere. The resulting fits for the four different 
types of crystals are shown as examples in Fig. 3. 

The final result for a set of runs under study yields 

Az = (0.063 ± 0.040) mm, (3) 
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thus assuring us of the very precise alignment of the MWPC's and the Csl 
Calorimeter. 



5 PIBETA detector active target 



Our regular active target is a cylindrical plastic scintillation counter with a 
length of 50 mm and a diameter of 40 mm. The counter is segmented into 
9 elements, as shown in a cross sectional drawing of the detector's central 
region in Fig. 4. A central target tube with a 3.5 mm radius is surrounded 
with four identical quadrant segments that compose an inner ring with 15 mm 
outer radius. The second, outer target ring is made out of four identically 
sized tubular segments that are rotated by 45° with respect to the inner ring. 
These 9 pieces are wrapped individually in aluminized Mylar foil making them 
optically insulated from each other. The parts are pressed together with a 
rubber band and the whole assembly is wrapped with a black plastic tape. 
Each target element is viewed by a miniature (8 mm diameter photo-cathode) 
Hamamatsu R7400U photomultiplier tube via a 60 mm long, tapered acrylic 
light-guide and is therefore acting as an independent counter. A photograph 
of the partially assembled target detector is presented in Fig. 5. 

The analog signal from each of the target segments is divided by custom-made 
passive splitters into two branches. One side is discriminated in the time- 
over-threshold discriminator (CAMAC PS7106); its outputs are digitized in 
a Time-to-Digital Converter (FASTBUS LCR1877 TDC), as well as counted 
in a scaler unit (CAMAC PS7132H) that is read out every 10 s. The second 
branch is connected to an Analog-to-Digital Converter (FASTBUS LCR1882F 
ADC), which is gated with a 100 ns event trigger gate. 

High voltages of the target segments' PMTs are adjusted using a 4.1 MeV 
7r + — > p^v^ stopping muon line as a reference until the 9 detector gains are 
matched when viewed on a digital oscilloscope. The final gain adjustment is 
done off-line using the ADC spectrum of stopping pions (see the next Section). 

The five central target detectors by design receive practically the entire stopped 
7r + beam, each counting at a comparable rate. The counting asymmetry (left- 
right or up-down) causes a shift of the centroid, that leads to a change in the 
integrated detector acceptance ratio. This is an easily measurable asymmetry 
even after the Poisson corrections for unresolved double hits and accidental 
coincidences with Michel decays. 
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6 Beam profile from active target ADC spectra 



The analysis presented in this Section is based on the PIBETA data collected 
in 50 successive production runs, corresponding to about 2 days of uninter- 
rupted data taking at the nominal ir + stopping rate of ~ 5 • 10 5 s _1 . 

Two sets of ADC histograms for each of our two beam triggers are defined 
for 9 individual target segments. These histograms are filled off-line with an 
additional requirement that no particle showers are detected in the Csl calo- 
rimeter, thus suppressing hadronic interaction events and delayed weak n + or 
ji + decays with a final state positron, either of which could affect the target 
energy depositions. 

Representative energy spectra for one target piece (No. 1) are shown in Fig. 6. 
The events in which the stopping pion (or thru-going positron) deposits its full 
energy in this target element are clearly distinguished in the high-energy peak. 
The energy spectra of all 9 target elements are empirically well represented 
by Gaussian distributions with square-root function tails: 



where pi are the 5 free parameters of the fit. The events with the stopping 
particle sharing energy between two or more target segments due to scattering 
and/or beam spreading fall into a low energy tail. The low energy tail is 
also populated by positrons from weak pion and muon decays which miss the 
calorimeter. The peak-to-tail amplitude ratio depends on the pion stopping 
rate, the ADC gate width, as well as on the divergence of the beam at the 
front face of the target. The comparison of the data with the target energy 
depositions calculated with GEANT3 Monte Carlo [8] suggests that the beam 
divergence p r /p z is small, ~ 1% (see also Sec. 8). 

We have extracted the total number of pions and positrons stopping in each 
target segment by fitting their high-end energy spectra with Gaussian func- 
tions, integrating the fitted peaks and normalizing the integrals to the Monte 
Carlo predictions. The results are summarized in Fig. 7. The top panel shows 
the target counting rates for the "7r-in-beam" trigger, while the bottom dia- 
gram corresponds to the "e-in-beam" events. The standard deviations shown 
are calculated from the uncertainties of the fitted Gaussian parameters and 
reflect both the event statistics and the quality of the fits. The central segment 
and the inner ring 7r + counting rates are determined with an accuracy better 
than 1%. The outer ring segments count the tails of the beam distribution 
with a 3-5 % uncertainty. 



f(E) = p x exp 



UE-P2 
A ?3 




+ 9(p 4 - p 5 E) - PsE, (4) 



9 



We have examined the resulting target counting rates by simulating the tar- 
get's response to the variety of lateral stopping beam distributions with the 
help of a simple Monte Carlo (MC) integration program. Good agreement 
between the Monte Carlo simulation and the measured rates is achieved by 
approximating the 3-dimensional beam profile with three separate Gaussian 
distributions in the three orthogonal variables x (beam left), y (up), and z 
(beam downstream): 



f(x, y, z) = iVo j 9(x - x) ■ exp -- 
9(x — xo) ■ exp 

i[y-yo 

exp 




(5) 



where the function 6{x) = for x > and 6{x) = 1 for x > 0, o x l and u x r are 
standard deviations of the horizontal x beam profile beam left (Bl) and beam 
right (Br), o y and a z are beam widths in the vertical (y) and axial (z) projec- 
tions, and the beam spot is centered at the point (xo, yo, zo). The asymmetry 
in the horizontal (x) beam profile is caused by a 4 mm thick carbon degrader 
inserted in the middle of the 7rEl beam-line. The momentum-analyzed pions 
and positrons have different energy losses in the carbon absorber, and are 
therefore spatially separated in the horizontal plane. The passive lead colli- 
mator in front of the PIBETA detector then blocks out most e + 's and yU + 's, 
minimizing the non-pionic beam contamination. 

All other attempted descriptions of the relative counting rates (e.g., a Lorentzian, 
piecewise parabolic functions, a box spectrum with or without exponential 
tails, etc.) resulted in far worse fits of the data. 

The five parameters of the f(x,y) distribution, namely <j x l, <j x r, <j y , x , and 
yo are calculated using the MINUIT minimization code [9] by requiring the 
best agreement between the distribution (2) and the data in Fig. 7 (mini- 
mum x 2 )- The SEEK and SIMPLEX search methods are used to suppress 
the dependence of the convergence on the uncertainty of the Monte Carlo in- 
tegration. SEEK minimizes the function using the Metropolis algorithm by 
selecting random values of the variable parameters, chosen uniformly over a 
hypercube centered at the best current value. SIMPLEX is a genuine multidi- 
mensional function minimization routine robust even to gross fluctuations in 
the function value. Both algorithms converged to the same solution. 

The best parameter values obtained in the fits are summarized in the second 
column of Table 1. The corresponding two-dimensional distribution is shown 
in Fig. 8, where it is superimposed on the sub-division of our target. The 
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agreement between the measured counting rates and the simulated rates is 
within one standard deviation, with the \ 2 P er degree of freedom of 0.9. 



7 7T+ beam profile from back-tracking tomography 



7. 1 Longitudinal n + stopping distribution 



The thicknesses of the beam-defining detectors, namely of the forward beam 
counter, the active degrader and the active target, were chosen to make the 
40.6 MeV incident 7r + beam particles stop close to the center of the target. 
Our GEANT3 Monte Carlo simulation of the longitudinal vertex distribution 
of decaying 7r + 's is presented in Fig. 9. The input to the Monte Carlo is the ir + 
momentum spectrum in Fig. I. A histogram of Monte Carlo z coordinates of 
7r + decay vertices is a Gaussian function with a width of o~ z = 1.69 ± 0.01 mm 
and a flat upstream tail integrating to 0.86 ± 0.05% events. The z position 
spread originates mainly from the energy straggling of stopping pions: the 
momentum spread of the incident beam contributes just 0.2 mm (or 12%) to 
the over-all axial distribution spread. The upstream tail represents n + decay- 
in-flight events and strong interactions with detector material. 



7.2 MWPC tracking and positional resolution 



The main decay channel of 7r + 's stopped in the target is the n + — > — > e + 
decay chain: 

7T+ -> fl+U^, (6) 

followed by 

H + -> e + v e v^. (7) 

The first reaction, namely the main two-particle pion decay, results in a mo- 
noenergetic fi + with a kinetic energy of 4.1 MeV that stops within the target 
detector because of its short 1.34 mm range. The muons subsequently decay, by 
producing so-called Michel positrons with a typical continuous /9-decay energy 
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spectrum characterized by its 52.5 MeV (= m M /2) end-point. A GEANT3 sim- 
ulation predicts that only ~ 0.48 % of the /^ + 's escape the target before coming 
to a complete stop. 

The suppressed n e 2 decay mode 

7T+ -> e + u e (8) 



proceeds with the branching ratio B We2 = 1.24 • 10~ 4 [10] and results in a 
monoenergetic 69.8 MeV positron that deposits an average of 4.5 MeV in the 
target. All average deposited energies, particle ranges and losses calculated 
in simulation assume a Xjy ~ 8 mm nominal beam spot size centered on the 
target. 

Michel positrons [Eq. (7)] and n — > ev positrons [Eq. (8)] are identified in 
the segmented plastic veto hodoscope PV via their characteristic energy de- 
positions and tracked in the pair of cylindrical MWPC's [11] that surround 
the target region. Total positron energies are measured in the Csl calorime- 
ter where they produce fully contained showers. The particles' energy spectra 
in the Csl calorimeter are used to discriminate between 7r + and /i + decay 
processes. The 7r e 2 positron vertices in the target therefore coincide with the 
7r + stopping points. The Monte Carlo root-mean-square (rms) of the Michel 
vertices turns out to be merely 0.04 mm larger than the 7r e2 rms. 

The directional resolution of the MWPC's has to be known before one can use 
back-tracking of the charged particles to find the real distribution of vertices 
in the target. Cosmic muon events are ideally suited for the resolution study 
of the MWPC's. Clean samples of cosmic muon events are collected during 
weekly beam-off cyclotron maintenance periods. The cosmic event trigger re- 
quires the coincidence of the cosmic muon veto counter and a Csl calorimeter 
signal above the "low discriminator threshold" (~5MeV). 

In the off-line data analysis two additional stringent software cuts are im- 
posed: (i) Csl calorimeter energy deposition > 200 MeV, and (ii) exactly two 
reconstructed hits in each MWPC. These constraints effectively remove any 
extraneous non-cosmic muon background. We reconstruct the cosmic muon 
track from a pair of hits in one MWPC and calculate the intersections of that 
track with the other chamber. Histogrammed differences between the calcu- 
lated x c , y c , z c and measured intersection coordinates x e , y e , and z e represent 
the positional resolutions of the chamber. The typical azimuthal angle reso- 
lution A0 and axial coordinate resolution Az for the MWPCi chamber are 
shown as an illustration in two panels of Fig. 10. The A<p rms is 0.74° and Az 
rms is 0.97 mm. 
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7.3 Back-tracking tomography: reconstruction algorithm 



The tomography algorithm relies on the back-tracking data obtained with 
the MWPC's. In the off-line analyzer computer program we define an 80 x 
80 x 80 mm 3 volume containing the cylindrical beam stopping target divided 
into elementary cubic cells of fixed size. We use initially 80x80x80 elementary 
cells, each with a volume of 1 mm 3 . For each reconstructed 7r e2 event a positron 
path length inside each elementary cell that is intersected by a MWPC track is 
found first. The analyzer code keeps a cumulative sum of tracks' path lengths 
for each individual cell. The high statistics data set ensures that the cubic cell 
containing the most n e2 vertices accumulates the largest path length sum. The 
distribution of path lengths therefore reflects the distribution of the decaying 
particles' vertices originating from the different volume elements [12]. 

The above- described algorithm determines accurately the center position (xo, yo, z ) 
of the x n + , y n + , and z n + distributions. While the shape of our n + stopping dis- 
tribution can be approximated by the product of Gaussian functions (Sec. 6), 
the real lateral distribution is more complicated, exhibiting outlying tails. 

The Monte Carlo (MC) method is then used to find the correspondence be- 
tween the input track vertex distribution in the target and the output path 
length distribution in the elementary cells. In the first stage one has to opti- 
mize the cell size. 

We generate the MC tracks from an asymmetric Gaussian input distribution 
(Eq. 6) uniformly into the solid angle covered by the MWPC geometry. The 
points at which the tracks intersect the cylindrical chambers are smeared next 
by the measured MWPC positional resolutions. We simulate the experimen- 
tally measured particle direction by finding the track through smeared inter- 
section points associated with the two MWPC's. Processing the large statis- 
tics of the MC tracks obtained in this way with the tomography algorithm 
described above we simulate the Xa>r, Syr an d ^ z t projections of the path 
length sums. The Gaussian fits to the simulated path length histograms are 
thus connected to the input distribution rms's. The dependence of the recon- 
structed width of Y,y T on the unit cell size is shown in Fig. 11. The input rms 
value of the vertical beam profile was 10 mm. From this study we deduce that 
for an optimal cell size of 0.3 mm a small correction of Aa y = 0.008±0.005 mm 
has to be applied to the extracted rms of the vertical (y) vertex distribution. 

The horizontal/lateral (x) and axial (z) rms values of the vertex profiles are 
treated analogously and the small correction factors (< 0.1 mm) corresponding 
to the chosen cell size are obtained. The two-dimensional "contour" plot in 
Fig. 12 shows the correspondence between the input vertex distribution in the 
horizontal coordinate x and the simulated path length histogram. 
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The MC tomography reconstruction is compared with the experimental data 
in three panels of Fig. 13. The agreement between the data points and the 
predicted projections is very good: by scaling the width of the input f(x,y) 
vertex distribution and studying the changes in the resulting chi-squares we 
estimate the uncertainty in the extracted cr x + and a y + to be ~ 0.05 mm. 
These results are summarized in the third column of Table 1. The details of 
the analysis are given in Ref. [13]. 



8 Pion beam profile: Dalitz event e + e vertices 



The selection criteria for 7r° Dalitz decay events, 7T° — > e+e 7, were: 

(1) no in-time hits in either active collimators or in cosmic muon vetoes; 

(2) no prompt 7r + in either the forward beam counter or in the active degrader; 

(3) total calorimeter energy deposition 10 MeV < E c < 200 MeV; 

(4) at least two minimum ionizing (MI) particle tracks pointing back to the 
target volume, identified as e+'s via energy depositions in the plastic veto 
and Csl calorimeter; 

(5) at least one neutral shower with E 1 > 20 MeV in the Csl calorimeter within 
±5 ns of two MI tracks; 

(6) if two or more candidate e ± and/or 7 tracks are found in a single event, the 
(7,e + ,e~) combination with the closest timing is selected. 

The candidate events satisfying the conditions listed above were written into 
a data summary tape in the off-line data replay. 

The e + and e~ tracks were back-projected to the target region and a pair 
of points, one on each track, which are closest to each other was found. The 
mid-point coordinates on the line connecting the closest two track points are 
defined as an event vertex. 

Fig. 14 shows the projection of Dalitz event vertices onto the horizontal (x), 
vertical (y), and beam (z) axes. The extracted peak widths, defined as a beam 
right and beam left Gaussian fits to the vertex distributions, are entered in 
Table 1, fourth column. 

The widths of the beam stopping distribution as well as its position are in 
excellent agreement with the tomography solution discussed in Sec. 7.3 above. 
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9 7r+ beam profile and 7r° — > 77 opening angle 



We have developed a complete GEANT3 Monte Carlo description of the PI- 
BETA detector that includes all major sensitive as well as passive detector 
components. The user code is written in a modular form in standard FOR- 
TRAN 77 and organized into over 300 subroutines and data files [14] . 

The lateral distribution of the 7T/3 event vertices in the x-y plane in our simu- 
lation is derived from the tomographic beam spot reconstruction of Fig. 13. In 
the direction of the beam axis the Gaussian n + stopping spectrum is generated 
with a z = 1.7 mm, as shown in Fig. 9. 

The MC 7r/3 event generator takes into account the phase space probability, 
the square of the matrix element and the radiative corrections [15,16]. The 
GEANT3 tt[3 detector trigger relies on the response of the Csl calorimeter; it 
requires two coincident neutral showers above the high energy discriminator 
threshold (HT) of 51.5 MeV that reconstruct into a tt° originating from the 
ir(3 decay at rest. 

The photon directions calculated by a logarithmic weighting method of shower 
localization [17]. The method is motivated by the exponential fall-off of the 
shower's transverse energy deposition [18]. The angular direction (6* c , <p c ) of 
a particle initiating an electromagnetic shower in the PIBETA calorimeter is 
found as an energy-weighted mean: 

N 

0c = ^ , (9) 

i=0 

where $i is the polar angle of an individual Csl module, ooi is a weighting func- 
tion, and iV is the number of nearest neighbors of the crystal with maximum 
energy deposition E . The sum is over the central crystal "0" and all of its 
nearest neighbors and thus contains 6 to 8 terms, depending on the shape of 
the centrally hit module. The formula for the azimuthal angle C is identical 
to that for the polar angle with the replacement 0, — > (pi. 

The weighting function Ui is defined by: 

Ui = max[0, a + ln(Ei) - \n(E tot )}, (10) 

where the best value for the constant ao is determined empirically to be 5.5, 
and E tot is the total energy deposited in iV + 1 modules that define the shower 
cluster. 
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The shape and absolute integral of the pion beta decay 7-7 opening angle 
histogram is very sensitive to the distribution of vertices of the initial decaying 
7r°'s. The acceptance A^J for the high discriminator threshold HT calculated 
in a GEANT3 Monte Carlo simulation is 

A*J = (1 - 0.4928 • KT 3 • r - 0.7365 • 10~ 3 • r 2 ) x 

(0.5977 - 0.2350 • 10~ 2 • s r . - 0.5500 • 10" 3 • s 2 r ), (11) 

where r = y/x 2 + y 2 is a lateral offset of the beam center, and s r is a lateral 
width (a standard deviation) of a Gaussian-shaped beam spot, where both 
variables are expressed in millimeters. This expression represents our main 
result, parameterizing the acceptance of the PIBETA detector. 

Fig. 15 shows 7-7 opening-angle histograms calculated using a GEANT3 code 
for three different beam spot sizes (top panel) and varying axial positions 
of the beam center (middle panel). The figure indicates that it pays to keep 
the beam spot root-mean square spread s r as small as possible, because the 
detector acceptance for a tight event vertex distribution is less sensitive to 
its shape and position. In particular, the integrated nf3 event acceptance is 
relatively insensitive to the longitudinal position z of the stopping point, 
with A^J changing <0.07% for z = ±10 mm, and is therefore neglected in 
Eq. 11. 

On the other hand, the resulting ir{3 acceptance is affected considerably by the 
spreading and lateral offsets of the beam profiles. Using the parameterization 
of Eq. (11) and propagating the uncertainties of the beam spot variables r 
and s r we obtain the results in Table 2. At our nominal stopped 7r + beam 
width of a X; y ~ 8 mm, the change in the rms value or lateral position of the 
beam spot by ~lmm affects the ir(3 acceptance by ~1.5%. Therefore, the 
experimental knowledge of the beam geometry at the sub-millimeter scale 
is required to accomplish the goals of the PIBETA experiment. This paper 
demonstrates that the required accuracy has been achieved. If absolute nor- 
malization is done with respect to the 7r e2 acceptance, the constraint on the 
beam geometry uncertainty is not changed significantly because in our exper- 
imental configuration A ne2 is almost independent of the small variations from 
the nominal beam model. The geometrical acceptance A ne2 varies less than 
0.5% for r , z , and s r in 0-10 mm range. 

In the bottom panel of Fig. 15 we present the MC # 77 calculation for our 
average beam distribution. The measured # 77 spectrum, based on an off-line 
replay of 59,011 7r(3 events is represented by full markers. For data taken over 
such a long time period we must take into account the temporal stability of 
the beam distribution. The back-tracing tomography reveals that our beam 
spot position and its size were remarkably stable. The reconstructed lateral 
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(rr) and longitudinal (zr) centroids of the stopping distribution vary with 
rms r ~ 0.04 mm and rms 2 ~ 0.2 mm, respectively, as shown in Fig. 16. The 
oscillations of the beam spot location and size are included in our GEANT3 
event generators. After that adjustment, nf3 data points and the MC simula- 
tion of # 77 agree with each other within the statistical uncertainties. 



10 Conclusion 



Precise knowledge of the PIBETA detector acceptance for the decay products 
of the stopped pion beam is required for the evaluation of the absolute branch- 
ing ratios of pion and muon decay processes. The geometrical acceptances of a 
compact detector like ours, defined with the specific cuts on the energies and 
angles of measured particles, depend strongly on the distribution of vertices 
of decaying particles. In this paper we show that we can measure position 
and the volume distribution of the stopped pion beam with sub-millimeter 
precision. The uncertainty in the 3-dimensional beam shape or position of the 
order of ~ 0.1 mm then translates into ~0.12% systematic uncertainty of the 
acceptance calculated with the GEANT3 Monte Carlo simulation. The con- 
sistency of our four different methods of extracting the pion beam stopping 
distribution is summarized in Table 1. 
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Fig. 1. The momentum spectrum of the beam 7r + 's entering the active de- 
grader /target counters calculated by the TURTLE program [6]. 

Fig. 2. A "7r-in-beam" trigger is defined as a four-fold coincidence between a forward 
beam counter (Ch. 1), an active degrader (Ch. 2), an active target (Ch. 3), and a 
19.75 ns rf cyclotron signal (Ch. 4). 

Fig. 3. Determination of the offsets between coordinate systems of the MWPC's 
and Csl calorimeter. The MWPC's track intersections with the calorimeter sphere 
for which four representative crystals also have the maximum energy depositions. 

Fig. 4. Cross section of the PIBETA detector central region (drawn in GEANT3). 
Shown are the segmented active target AT, a pair of cylindrical MWPC's, and the 
plastic veto hodoscope array PV. 

Fig. 5. A photograph of the segmented PIBETA active target. Tapered acrylic 
light-guides are glued to the mutually optically isolated target segments. 

Fig. 6. Representative energy spectra of beam pions (7r-in-beam trigger, top panel) 
and beam positrons (e-in-beam trigger, bottom panel) stopping in the target seg- 
ment number 1. Monte Carlo GEANT3 predictions are represented by full line 
spectra. 

Fig. 7. Integrated counting rates in 9 target segments for the 7r-in-beam and 
e-in-beam triggers, respectively. Orientation of the coordinate system and target 
segment labels (0-8) are also indicated. 

Fig. 8. The 2-dimensional shape of the 7r + beam superimposed on an outline of the 
segmented target. The reconstruction is based on the target ADC spectra collected 
during years 1999-2001. 

Fig. 9. The z-component of the tt + stopping distribution in the active target calcu- 
lated in the GEANT3 simulation using the momentum distribution in Fig. 1. 

Fig. 10. The MWPCi directional resolutions in the azimuthal angle (f> (top panel), 
and axial coordinate z (bottom panel) extracted from a sample of cosmic muon 
events. The data points for <f> and z coordinates are fitted with a sum of two Gaus- 
sians, the broader one accounting for the cosmic muon scattering in the apparatus. 

Fig. 11. The dependence of reconstructed vertical rms widths a y of the beam stop- 
ping distribution on the size of the elementary cell used in the calculation. 

Fig. 12. A lookup table for determining (T xL j R (L left, R right) is calculated starting 
from the synthetic source distributions in the Monte Carlo. The fitted "contour" 
lines correspond to the values of equal cr's. 
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Fig. 13. The distribution of path length sums for x (top), y (middle) and z (bot- 
tom) coordinate projections. The experimental points for back-tracked ir + — >• e + is e 
positrons are compared with the Monte Carlo simulation that uses the beam stop- 
ping profile deduced from back-tracking tomography. 

Fig. 14. The horizontal (x, the top panel) and vertical (y, the middle panel) beam 
profiles deduced from the Dalitz event 7r° — > je + e~ vertices. The axial z coordinate 
cut around the target center, Az = ±10 mm is imposed to select the tt + stopping 
distribution (the bottom panel). The agreement with the tomography solutions is 
very good. 

Fig. 15. The 7-7 opening angle following the 7r° decay from the tt + — ► 7r°e + f process 
as reconstructed in the PIBETA segmented calorimeter. The GEANT3 prediction 
for the stopping n + beam profiles with three different transverse widths (top panel) 
and three different longitudinal stopping points (middle) is shown. The bottom 
panel shows a comparison of the partial nf3 events data set with the Monte Carlo 
prediction. The plotted uncertainties are statistical only. 

Fig. 16. The spatial stability of the n + beam spot center during one three month data 
taking period. The lateral beam center location (top) and the axial beam stopping 
point (bottom) are reconstructed via back-tracking tomography. The reconstruction 
uncertainty for a single run is <0.02mm. 
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Table 1 

The rms widths and centering of the 7r + beam spot in the PIBETA active stopping 
target deduced by four different methods: (i) using the counting rate matrix in 
the segmented target, (ii) via back-tracking of tt + — ► /i + — ► e + positrons with two 
cylindrical MWPC's, (iii) from Dalitz decay e + e~ pair vertices, and (iv) from fits 
to the opening angle distribution of two tt° — ► 77 photons. BR/BR stands for beam 
right/left, BU/BD labels beam up/down directions. 



Beam profile 


Segmented 


MWPC's 


Dalitz e + e 


— ^'yy 


parameters 


target rates 


back-tracking 


vertices 


opening angle 



(mm, BR horizontal) 
(mm, BL horizontal) 

<7y (mm, BU vertical) 

a® (mm, BD vertical) 

a z (mm, axial) 

x (mm) 

Do (mm) 

z (mm) 



7.0 ±0.3 
8.8 ±0.3 
7.0 ±0.3 
7.0 ±0.3 

-1.0 ±0.3 
0.1 ±0.3 



7.04 ±0.04 
9.50 ±0.07 
7.44 ± 0.06 
7.46 ± 0.07 
1.70 ±0.03 
-1.92 ±0.05 
-0.12 ±0.05 
0.30 ±0.02 



6.67 ±0.37 
8.46 ± 0.75 
7.25 ± 0.49 
7.40 ± 0.32 

1.8 ±0.3 
-0.95 ±0.70 
-0.31 ±0.60 
0.43 ±0.12 



7.0 ±0.5 
9.0 ±0.5 
7.6 ±0.5 
7.6 ±0.5 
2.0 ±0.3 

-1.0 ±0.3 
0.0 ±0.3 

0.30 ±0.1 



21 



Table 2 

Summary of two main uncertainties connected with the size and position of a ir + 
beam spot in determining the detector acceptance A^p. Relevant geometrical vari- 
ables are a lateral offset of the 7r + beam center tq with respect to the geometrical 
center of the target and a lateral rms width of the ir + beam profile s r . Beam spot 
parameters obtained from the back-tracking tomography are used. 





Source 


General 
value (%) 


Current 
value (%) 


W)- 1 


' (9A^J/dr ) ■ a ro \ ro=lmm 
■ {dA^J/ds r ) ■ CT Sr | Sr=8mm 


1.15 ■ a ro (mm) 
1.26 ■ a Sr (mm) 


0.06 
0.10 
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